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Abstract 

The dispersion relation in matter of ultrahigh-energy neutrinos above the 
pole of the W resonance (E v > 10 7 GeV), is studied. We perform our calcu- 
lation using the real-time formulation of Thermal Field Theory in which the 
massless limit for the W boson is taken. The range of active-to-sterile neu- 
trino oscillation parameters for which there is significant mixing enhancement 
during propagation through the interior of the Earth, and therefore significant 
attenuation of neutrino beams in the Earth at high energies, is estimated. Fi- 
nally, this range is considered in view of the cosmological and astrophysical 
constraints. 
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It is now well established |]-[3| that the Earth's diameter exceeds the attenuation length of 
neutrinos with energies greater than 25 TeV. Such an estimate was based on the calculation 
of the cross sections for vN collisions at ultrahigh energies (UHE), [E v > 1 TeV). Because 
of the smallness of the electron mass, ve interactions are generally considered as negligible 
with respect to vN interactions; vN interactions therefore provide the dominant signal and 
account for most of the attenuation of neutrino beams in the interior of the Earth at ultrahigh 
energies JT^. There is one exception though, the resonant formation of the intermediate 
W~ boson in ve interactions in the neighborhood of E r v es = My V /2m e ~ 6.3 x 10 15 eV. 

The promising tool for detection of UHE cosmic neutrinos by means of neutrino telescopes 

consists of recording the long-range muons produced in charged- current vN interactions 
that occurs in matter surrounding the detector. Apart from efficient shielding from the flux 
of atmospheric muons, such upward-going muon events have the advantage of enhancing the 
effective volume in proportion to the range of the produced muons (typically a few kilometers 
for ~ 10 TeV). For our purpose it is important to note that the rate for upward- 
going muons does not depend only on the probability for neutrino conversion to a muon 
with energy above the threshold energy, but also through the interaction length which is 
responsible for the attenuation of the neutrino flux due to interactions in the Earth's interior. 
The typical situation that occurs for E v < 10 5 GeV is that the upward rates depend little 
[ ij on the calculated vN cross sections, since the enhanced (weakened) interaction rate is 
nearly compensated by the enhanced (weakened) attenuation of UHE neutrinos propagating 
through the Earth. On the other hand, the detection of cosmic neutrinos at energies of 
10 16 eV or larger is beset by the problem of the increased importance of attenuation of 
neutrino beams [B10. However, even so, the upward rates produced by neutrinos from 
powerful radiation sources, like Active Galactic Nuclei (AGN) ||, should be observable in a 
detector whose effective area is A ~ 0.1 km 2 . 

As for event rates involving electron neutrinos, they are generally smaller than the muon 
event rate by the flux ratio (the initial fluxes of UHE neutrinos originating from AGNs are 
expected to have a ratio v e jv il ~ 1/2) times the detector length divided by the mean muon 
range, because of the rapid energy loss of electrons (or annihilation for positrons). Still, it was 
shown recently that the Landau-Pomeranchuk-Migdal effect ]7[ may effectively enhance 
the electron range by detecting upward-going air showers initiated by the v e interaction 
near the Earth's surface. On the other hand, resonant ve scattering contributes significantly 
to the attenuation of z/ e 's, meaning that the flux of electron antineutrinos in the range 
2 x 10 15 eV < E u < 2 x 10 16 eV is extinguished for neutrinos traversing the Earth 

In the present paper we are going to consider another mechanism for attenuation of UHE 
neutrinos propagating through the Earth, namely, matter enhanced neutrino oscillations 
^ e ^ Vs, where s is a sterile neutrino (e.g. a singlet under the gauge symmetry of the 
Standard model). We shall be concerned exclusively with the case of v e <-> v s oscillations 
where the energy of UHE neutrinos is above the resonant energy, E u > E r u es . Owing to 
the new form of the effective matter potential in this regime, we are in position to study 
MSW resonance effect || previously ignored in the literature. On the other hand, the 
matter effect of the Earth through the standard effective potentials [E v << El es ), in the 
region of oscillation parameters relevant for the solar and atmospheric neutrinos as well as 
those suggested by the LSND result, is well established now ||. Even more, a new effect 
of matter-enhanced neutrino mixing, based on a maximal constructive interference among 
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transition amplitudes, has been discovered recently [flO | . In the following, we shall first 



derive the induced mass squared of the electron neutrino with E u > E r v e& then we find the 
range of neutrino parameters for which there is matter-enhanced v e «-> v s oscillation during 
propagation of v e through the Earth and finally we discuss if the established range could 
survive constraints from type II supernovae as well as big bang nucleosynthesis on v e <-> v s 
mixing. 

Effects of a medium on neutrino propagation is determined by the difference of potentials, 
whose standard-model contribution in the context of Thermal Field Theory (TFT) may 
readily be obtained from the relevant thermal self-energies of a neutrino: charged- current, 
neutral-current and tadpole. Of these, only the charged- current and tadpole contribution are 
relevant for our consideration. Let us consider the charged-current diagram in more detail. 
Since the "target" electrons are massive, the W boson may be considered "massless" always 
when s ~ 2E u m e 3> M^. This corresponds to energies in the lab frame E u >Qx 10 15 eV. 
In the opposite limit, s <C M^, the W boson should be considered "massive" and the usual 
contact approximation for the ^-propagator is adequate. Using the real-time formulation 
of TFT, we discover by explicit calculation that the induced mass squared of u e is equivalent 
to the fermion thermal mass squared JTT| , 
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A » = 2^ I- kdkn ^ iE,.»E;r). (i) 
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where g ~ 0.63 is the gauge coupling constant. Note that in contrast to the mass squared 
induced by the standard MSW potential, 

A c : = 2V2G F N e E u [E v < E r u es ) , (2) 

(1) is independent of neutrino energy and also there is no explicit dependence on the number 
density of electrons. It should be clearly stated here that actually we are not dealing with 
field theory in equilibrium since all the electrons in the medium are bound electrons. Still, 
one is allowed to retain the usual real-time formalism by taking the (ll)-component of the 
electron propagator to be 



S 11 (k) = (}6 + m e )(- 2 1 - T - - + 2mn e (k Q )5(k 2 -ml)) , (3) 



where now a bound electron is assigned a distribution n e (k ). Thus (1) describes the plane- 
wave impulse approximation, which is, for instance, the basic approximation of electron 
momentum spectroscopy of atoms and molecules [12|. 



Although a connection with equilibrium TFT is now established, one may still wonder 
why (1) is equivalent to the fermion thermal mass, that is, to the characteristic mass scale 
which appears naturally only in the high-temperature limit of the fermion self-energy. Before 
going into details, let us stress that to first order in perturbation theory at high temperature, 



the poles of the full fermion propagator are determined just by the thermal mass [IT]. In 
hot theories, the high-temperature limit means that the temperature is much larger than 
mass of the particles under consideration (m e , My/ in our case) and the external momenta. 
Consequently, all particle masses can be ignored for practical purposes. In the standard 
electro-weak theory we now show that the opposite case, when the external momenta are 
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large and the characteristic scale of the medium is much less than the particle masses, can 
faithfully mimic the high-temperature limit. Indeed, when E v 3> E r v es , the four-momentum 
squared of the W boson is always much larger than Mj^, and hence the W boson can be 
considered "massless" . To get rid of the electron mass, note that it appears in the numerator 
of the thermal part of the electron propagator, and both in the numerator and denominator 
of the vacuum part of Sn [see Eq. (3)]. It is however trivially to see that m e from the 
numerator disappears when sandwiched between the two electroweak vertices 7 M L, where 
L — ^(1 — 75) is the projection operator for the left-chiral fermions. On the other hand, the 
real part of the neutrino self-energy is given only by a contribution where a cut is through 
the electron line, since W bosons are absent from the medium. This means that there is 
no contribution from the vacuum part of Sn, where m e appears in the denominator . This 
completes the proof that the high-energy limit is equivalent to the high-temperature limit 
of hot gauge theories. Strictly speaking, since the thermal mass always involves the thermal 
effects from both the W boson and the electron propagators, the only difference is that the 
former is absent in (1). 

Using the above interpretation for bound states and keeping the same normalization as 
for quasifree states, one can rewrite (1) in the following form, 

A? ~ 0.2 < k' 1 > N e , (4) 

where < k > is the average momentum of bound electrons. For the rough estimates presented 
here, it is sufficient to assume that < k^ 1 > ~ < k > -1 . Let us choose the average 
momentum per atom (with the atomic number Z) as a quantity of interest here. Going 
back to atomic physics, one can determine this quantity by applying the Thomas-Fermi 



method [13] to the calculation of the total ionization energy of a neutral atom. The result 



is 



< k > z ~ 4.6 Z 2/3 keV . (5) 

For our purpose, let us recall that the interior of the Earth consists of two regions of slowly 
varying density - the core and the mantle, with particularly strong density change between 



the lower mantle and outer core. The density profile of the Earth can be found in Jm 
The density of the mantle increases from 3 to 5.5 gcm~ 3 (the average value is 4.7 gem -3 
and the average electron fraction is 0.49), while the density of the core varies from 10 to 
13 gem -3 (the average value is 11.8 gem -3 and the average electron fraction is 0.47). The 
core comprises heavier elements, presumably nickel (Z = 28). Hence from (5) we have 

<k>l re ^42keV . (6) 

The mantle consists of lighter elements [Z = 8 — 16), and our estimate in this case is 

< k > mantle — 25 keV . (7) 

Before determination from the resonance condition of a range of neutrino masses where 
maximum mixing enhancement may occur, one should consider the tadpole graph as well. 
Being a constant independent of the external neutrino momentum, it is the same as in the 
standard MSW case [E v <C ££ es ), giving rise to the induced mass squared which grows 
linearly with energy, i.e., 
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A tadpole = _^2 GpNnEi 



(8) 



where N n is the neutron number density. Apart from a negative sign in (8), let us com- 
pare the magnitude of (8) with the charged- current contribution (1). It turns out that 
for E v > 10 9 GeV (8) is beginning to dominate over (1). For our purpose, it is however 
enough to consider E v ~ 10 7 — 10 8 GeV since the short v e interaction length for energies 
E v > 10 9 GeV means that the flux of electron neutrinos is extinguished for neutrinos travers- 
ing the Earth. Hence for E v < 10 9 GeV, the total induced mass squared for v e is essentially 
given by (1). It is interesting to note that at a particular energy around 10 9 GeV there is a 
nearly complete cancellation of matter effects in the neutrino propagator due to the sign of 
(8). ' 

With the above simplifications and the range of densities in the Earth's interior as dis- 
cussed before, one finds from the resonance condition, A c v c pa Am^, a range of neutrino 
masses where maximal mixing enhancement may occur, 

0.07 < AmlJkeV 2 < 0.12 (10 9 GeV > E v > 10 7 GeV ). (9) 

Here we have taken a small mixing angle, cos 29 ~ 1, in order to study oscillation en- 
hancement. Am^ s > in (9) means that v s is heavier than v e . This is also true for v e <-> v s 
oscillations as A u c > (in contrast to the standard MSW potential where the sign is reversed 
for antineutrinos). 

Notice that because of the resonance condition which is energy independent and the 
range of densities in the Earth, the range (9) is actually very small. Even so, it lies in the 
region which might be very interesting to astrophysics as well as cosmology. We recall that 
the prediction for the sterile neutrino of ~ keV mass is not in contradiction with any of the 
present bound. Indeed, the ~ keV mass is needed if active-to-sterile neutrino oscillations 
are to solve the pulsar velocity puzzle fL5f . In contrast to active-to-active oscillations, this 



solution is not in conflict with the cosmological bound on stable neutrino masses since the 



~ keV mass sterile neutrino has been proposed as a viable dark-matter candidate [|16[ . 

It is well known that if we add the reported results from the LSND collaboration [17| to 
the list of neutrino anomalies, then the explanation of all of them requires a four-neutrino 
scheme with three active neutrinos u e , u^, v T and one electroweak-singlet neutrino ||18|| . This 



introduces the 'LSND gap' of a few eV, a few orders of magnitude below our preferred range 
(9). However, the LSND result is not confirmed (although not completely ruled out) by a 
similar KARMEN experiment |19[|. 



It is easy to estimate the range of active-to-sterile neutrino parameters for which there 
is significant enhancement mixing during propagation through the Earth. For significant 
transitions to developed, it is necessary that the propagation distance be greater than about 



a quarter of a wavelength at resonance p0| . This constraint gives us a lower limit on the 

trough the Earth (2R earth ) we have, 

< ZRearth , (10) 



mixing angle. Taking the longest distance through the Earth (2R eart h) we have 

irE„ 



Amg S sin 26 

which for E v = 10 7 GeV and 2R earth = 1.27 x 10 9 cm gives (sin 26) core > 2 x 10" 3 and 
(sin20) mantte > 3.4 x 10~ 3 . 
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There is however a stronger limit on sin 29 coming from the condition for unsuppressed 
oscillations. The oscillation frequency must be real, otherwise the system is critically over- 
damped, the oscillations would be fully incoherent, and hence in fact there will be no oscilla- 
tions. Since for neutrino energy in the range 10 15 eV < E v < 10 21 eV, the cross section scales 
with E v as a oc E® A one finds the mean free path for neutrinos with E v = 10 7 GeV to 
be I — 0.1 x 2R eart h- The condition for unsuppressed oscillations, 

l m <2l, (11) 

then gives (sin26') core > 4 x 1CT 2 and (sin 26) mant i e > 7 x 1CT 2 . Notice that although the 
Earth is opaque to UHE neutrinos, the oscillations may proceed unsuppressed whenever the 
above requirement is satisfied. 

Let us finally check up the range (Am^/fceV 2 ~ 0.12, sin26 l > 4x 10~ 2 ) from a viewpoint 
of astrophysics and cosmology. The effect of a resonant v e v s mixing on a type II 
supernovae was considered in [21]. The bounds on (Am 2 s ,sin2#) derived in [21] are valid 



only if the sterile neutrinos have a mean free path larger than the radius of the supernova 
core after passing the resonance; this is the case if sin 29 < 3 x 10~ 2 . Our range is therefore 
unaffected by the type II supernovae constraint. On the other hand, a naive bound on the 
v e «-> v s mixing from big bang nucleosynthesis was derived, Am 2 s sin 4 2# < 5 x 10 -6 eV 2 



22j • Notice a disagreement of our preferred range with the above naive bound. However, 
the naive calculations ignored the creation of v — v asymmetries by active-sterile oscillations 
23fl in the early universe; these may efficiently suppress v e <-> v s oscillations, and therefore 



invalidate the conclusions drawn from naive calculations (even maximal <-» v s oscillations 



as a solution of the atmospheric neutrino anomaly p4j cannot be excluded |25|). The 
measurement of the flux of UHE neutrinos could thus provide us with a new test for this 
cosmological scenario. 

To summarize, the dispersion relation for electron (anti) neutrinos in the Earth's interior 
for energies above the pole of the W resonance, is derived. Then we have considered MSW 
oscillations for cosmic neutrinos traversing the Earth by including the charged-current self 
energy diagram for u e . We have shown that the range of neutrino masses where maximal 
enhancement may occur could be interesting from a viewpoint of astrophysics and cosmology. 
Let us finally stress that in order to study a nadir angle dependence beyond 34°, where 
neutrinos always propagate outside the core, a weaker attenuation of a v e beam would require 
the inclusion of the tadpole self-energy for energies beyond 10 9 GeV. This interesting case 
in now under study. 
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